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Abstract 

We describe a new method in lattice field theory to compute 
observables at various values of the parameters Aj in the action 
S[4>, Aj]. Firstly one performs a single simulation of a "refer- 
ence action" 5*[0 r ,A[] with fixed A[. Then the ^-configurations 
are transformed into those of a field distributed according to 
S[(j), Aj], apart from a "remainder action" which enters as a 
weight. In this way we measure the observables at values of Aj dif- 
ferent from A[ . We study the performance of the algorithm in the 
case of the simplest renormalizable model, namely the 4 scalar 
theory on a four dimensional lattice and compare the method 
with the "histogram" technique of which it is a generalization. 
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1 Introduction 



It is often necessary to study the observables of a given lattice field the- 
ory over a wide range of values of the parameters (coupling constants and 
masses). This is the case, for instance, when mapping out the phase structure 
of the model [0 or in problems which require a fine tuning of the parame- 
ters 0, or in the study of the lattice gauge theory beta-function ||. In 
a numerical approach to these problems one should perform independent 
Monte Carlo simulations for each desired value of the parameters. If there 
are many couplings and masses this method may become prohibitively time 
consuming. 

The so called "histogram" technique [|J is usually adopted to deal with 
these problems. Suppose that one wants to study observables for a lattice 
theory of a field with action S[<f>, AJ for various values of the parameters Aj. 
One runs a single Monte Carlo simulation using a reference action S[4>, A£] 
with a fixed set of values of the parameters A[. Then, by reweighting the 
generated configurations with the factor e AS , where AS = S[<f), A[] — S[(p, AJ, 
one computes the observables associated to the action S[(f>, A*]. 

This method is subject to various limitations ||, H|. First of all the fluc- 
tuations of AS* are of the order of the square root of the lattice volume. 
Therefore severe statistics problems are met, unless the lattice volumes are 
not too large. More importantly, the statistical significance and the relia- 
bility of the results depend on the size of the overlap of the regions of the 
configuration space covered by the importance sampling procedure based on 
the actions S[4>, A£] and S[4>, A»] respectively. As a consequence, the reliability 
of the results is not known a priori [^, ||. 

We propose a new method to deal with these problems. Suppose that 
one finds a "field transformation" from 0, distributed according to the action 
S[(p,Xi), to a "reference field" r , distributed according to the "reference 
action" S[(j) r , A£]. More precisely the fields and r are related by 

D0 exp{-S[0,A 4 ]} ~ D(f) r exp{-5[^,A[]}, (1) 

with Dip the usual (naive) field measure. Then by a single Monte Carlo 
simulation with action S[(p r , A[] one generates (^-configurations which can 
be transformed into the importance sampling configurations for 0. 
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In general, it is difficult to construct exactly the transformation (|1]) and 
one has to resort to approximations. In this case one is left with a "remain- 
der action" 5S[cj) r ] which enters as a factor e ss in the right hand side of (|I|). 
Therefore, also in this approach, one needs to reweight the (/^-configurations 
generated in the Monte Carlo run. The efficiency of the method depends then 
on: i) whether a simple enough field transformation can be constructed; ii) 
how large are the fluctuations introduced by the reweighting from the re- 
mainder action, iii) how large is the overlap of the 0- and of the trans- 
formed (//"-configurations. The field transformation method is a generaliza- 
tion of the "histogram" technique which, in our language, corresponds to 
choose the trivial mapping r = 0, so that the remainder action is given 
by A,S[0] = S[<j>, \j] — S[4>,\i\. The advantage of our method is that one 
can devise a systematic procedure for improving the field transformation by 
reducing the "remainder action" and therefore by increasing the size of the 
overlap of the <fi and of the transformed <p r configurations. 

We shall study the efficiency of the method in the case of the the simplest 
renormalizable field theory: the 4 model on a four dimensional lattice with 
L 4 sites. The action depends on two parameters, A 2 and A 4 , the coefficients of 
the quadratic and quartic field monomials respectively. As a first example we 
shall consider an ultralocal field transformation, i.e. one which is the same 
for all lattice sites. We shall study the magnetic susceptibility x an d the 
second moment of the correlation function /z 2 , for various values of the lattice 
parameters near the critical line where the square mass gap m 2 = 8x//J>2 
vanishes. We shall consider lattice sizes L = 8, L = 16 and L = 20. 

Particular attention will be devoted to the statistical fluctuations coming 
from the reweighting due to the remainder action. Since the field transfor- 
mation method is a generalization of the histogram method, the analysis of 
the errors can be done by following the same procedure as in Ref. |5|, 

Comparing the performance of the two methods, we find that, for the 
region of parameters and lattice sizes considered, the field transformation 
method always gives accurate values for the observables. The histogram 
method gives estimates of comparable accuracy only when the quartic cou- 
plings are the same, A4 = A 4 . When this is not the case, the estimated values 
of the observables are not accurate although the remainder action has no 
large fluctuations. 
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In Sect. 2 we describe the method for the scalar model. It will be clear 
that the method can be generalized to other models. In Sect. 3 we describe 
the numerical performance of the simplest field transformation. In Sect. 4 
we comment our results. 



2 Description of the field transformation 
method 

We consider the scalar field theory on a four dimensional lattice with periodic 
boundary conditions and L 4 sites. The lattice action is 

S[(j), A 2 , A 4 ] = (fin+fi + 4>n-n) + V ( ( f ) n) 

nfi n 

= ^V*+£#. (2) 

Only two parameters A 2 and A4 are needed. The normalization of the field 
is fixed by the kinetic term. The field 4> n is defined on the lattice site at 
the position given by the four-vector n and the sum over fi extends to four 
dimensions, \i — 1, • • • 4. 

The Green functions, given by the expectation of field polynomials 

P[0] = <f>m <Pn 2 ■■■(f>n k , 

are defined by 

(P[0]> = I / J] d^ n e~ s ^ P[0] , Z = /II #n e" 5 ^ . (3) 

Z J n J n 

Our aim is to evaluate the Green functions on a range of values of the lattice 
parameters Aj by using the information provided by a single Monte Carlo 
simulation with the reference action S[4> r , X^, A£]. To this end we look for 
a mapping of the field n onto a reference field 0^ such that (0) is at least 
approximately satisfied. The transformation is defined by the matrix 

j dj>n_ () 
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The field <p n is a functional of the reference field, <p n = n [0 r ], so that also 
the monomial P[<f>] is a functional of the reference field 

vw\ = KW'-'KW = P[^11- (5) 

By changing the fields according to (|4]), one has 

IJ#n e- s ^ = l[df n e~ s ^ e TrlnJ = l[df n e~ s ^ e 5s ^ , (6) 

n n n 

where the "remainder action" 5S, given by 

5S[(j) r ] = -S[<f>,\] + S[(j) r ,\l] + TrlnJ (7) 

allows for the fact that the transformation defined in (§) may not exactly 
satisfy (0). 

The Green function (|3|) can be expressed in terms of weighted expecta- 
tions as follows 

(VW] e 5s ^Y 
W]> = V [jsmy ' . (8) 

where (■ ■ -) r is the expectation value computed with the action S[<p r , A^]. 

To compute the Green functions for various values of A2 and A4, one 
proceeds as follows: by performing a single Monte Carlo simulation with 
S'[0 r ,A[], one generates a (sufficiently large) sequence of (^-configurations. 
Then, for each desired value of A2 and A4, one computes, by using (f|) and 
(0), the values of V[4> r ] and 5S"[0 r ] on each configuration. Averaging over 
the configurations one obtains the expectation values in the numerator and 
denominator of (§) and therefore the Green function. In the following we 
discuss the form of the matrix J in (|J) that we have considered. 

2.1 Form of the field transformation 

As stated in the introduction, an appropriate field transformation should 
satisfy two conflicting requirements: i) it should be easily computable; ii) it 
should produce a weight factor e 5S with small fluctuations. 
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Let us consider, as a first and simplest example, the following ultralocal 
field transformation 

_ d(j) n _ v{(t>n) _ v r {(f> r n) 

where C is a constant and t> and v r are given in (E3) with parameters Aj and 
A£ respectively. Implementing the initial condition n = for 0^ = the 
solution of (0) is 

f n d<\> = C [ K d<j> e~ vr ^ . (10) 

Jo jo 

Since both fields are non-compact we require that (p r n — > oo as n — >■ oo. This 
fixes the constant C as the ratio of the two integrals extended to infinity and 
defines <h n strictly monotonic function of (jf n . 

From (P) we find (neglecting a field independent term) 

TrlnJ = £K0 n )-^)], (11) 

n 

which gives (see (0)) as remainder action the difference of the kinetic terms 

53 = jE(W 0n +M + 0n-„) - ft( + &_ M ) ) • (12) 

This is the logarithm of the weight to be used in (^) for computing the Green 
functions. It is clear that too large fluctuations of the weight e ss would spoil 
the efficiency of the method. 

In order to analyze SS we consider its expansion around small values of 
the fields. From flip]) we obtain the following expansion of <ft n 

K = C<Pl (1 + O((0;) 2 )), (13) 

which gives 

SS = \{C 2 - 1) £ (&( ^ + + O ((« 4 ) (14) 

One expects that the first term quadratic in the fields gives the largest con- 
tribution to the fluctuations. Then one may try to improve the field trans- 
formation in order to be left with a SS not containing quadratic terms in the 
fields. 
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As already observed, our method reduces to the usual histogram 
reweighting technique if one chooses the trivial mapping (ff n = <p n . The cor- 
responding remainder action becomes the difference of the potential terms 
8S = £ n [ w (0n, Ai) - v((f) n , XI) }. 

3 First application 

Our aim in this section is to test how effective is the field transformation 
method by reproducing some standard results for the four- dimensional lattice 
model of the previous section. 

By using the simple field transformation defined in @ we compute for 
various values of Aj the following two physical observables: the susceptibility 

x(A 2 ,A 4 ) = jjY,(<Pn<Pm) , (15) 

*-* n,m 

and the second moment of the correlation function defined by 

/^(A 2 ,A 4 ) = ( n ~ m ) 2 (&» <Pm) • (16) 

^ n,m 

In terms of these quantities, the mass gap squared is defined by m 2 = 8x/l^2- 

There are various points to be explored in order to test the efficiency and 
accuracy of the field transformation method. The remainder action is the 
crucial quantity to be kept under control. Indeed too large fluctuations of 
5S will ruin the efficiency and the accuracy of the method. The issues we 
consider are the following: 

• the main expected challenge is the increase of the number of degrees of 
freedom, namely we might ask whether the quality of the results will 
worsen going to relatively large lattice volumes. The remainder action 
SS in ( P^D being a sum over the whole lattice, its fluctuations might 
grow rapidly with the square root of the volume and the method might 
become inefficient for large lattices; 

• at a fixed lattice size and for a given reference point A[ we should 
analyze as many values of Aj as possible, even relatively far away from 
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A[. Since the (absolute value of the) remainder action SS increases 
with the differences | A2 — A2 1 and | A4 — A4 1 (see (|HD), the efficiency and 
accuracy of the method might worsen dramatically by moving too far 
away from the reference parameters A^ - ; 

• finally we should test the method for Aj close to the critical line 
(A™ 1 *, A™ 1 ) of vanishing mass gap. Here one expects that the remain- 
der action would exhibit the large fluctuations which are typical of any 
quantity near a critical point. 

Most computations have been performed on a 8 4 lattice, which is suf- 
ficiently small to make a fast exploration of many values of Aj possible. A 
few computations on 16 4 and 20 4 lattices have been performed in order to 
explore the performance of the method as a function of the lattice volume. 

The region of parameters explored in our analysis for the 8 4 , 16 4 and 
20 4 lattices is represented in Fig. [l]. The reference point r = (A 2 , XQ = 
(—0.18,0.5), represented as a star, indicates the values of parameters A 2 , A4 
used for the reference action. The points in the regions denoted by A, B 
and C correspond to the values of parameters at which we have obtained 
the observables by using the field transformation method. At some of these 
points (black circles) we also have performed independent Monte Carlo test 
runs. The continuous line represents the one loop estimate of the critical line 
(A| rit , A™ 1 ) which, in the thermodynamic limit, is associated to the vanishing 
of the mass gap and to the divergence of the susceptibility. When approaching 
the critical points one has (L^oo) 

MA 2 , A 4 ) ~ x 2 (A 2 , A 4 ) -> 00 , m 2 = 8^- - , (17) 
3.1 Field transformation 

A good accuracy is necessary in the calculation of both the constant C and 
the functional dependence of 4> n on <p r n defined by fllPf). This is obtained by 
computing the two integrals F(<f>) = dt e~ v{t) and F r (0 r ) = jf dt e"^. 
The constant C is obtained by C = F(oo)/F r (oo) and is a function of Aj 
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Figure 1: Parameter space. The continuous line indicates the one-loop crit- 
ical line A™ 1 * = —0.46479 A| nt . The point r denotes the parameters (\ r 2 , \ r A ) 
of the reference action used in the Monte Carlo run. The performance of the 
field transformation method is studied by evaluating the observables for sets 
of points on the segments A, B and C (open and black circles). The black 
circles represent points at which Monte Carlo test runs have been performed. 
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and \\. For small difference of parameters it has the form 

(A 2 + 8) 2 (A£+8g \ vry^ 1 

y% >AJ )v(l/A) 2 + "\ ^ 

Then for any value of r , the equation F{$) = CF r ((p r ) is readily solved 
for and a table of correspondence of values for the mapping cff — > is 
created. This table is then used every time a measurement is performed 
for the mapped configuration 0. The increase in CPU time is negligible 
compared with the direct Monte Carlo simulation. 




3.2 Monte Carlo simulation 

We have performed several Monte Carlo simulations with the reference ac- 
tion, corresponding to the point r in Fig. |], on lattices of various sizes. We 
have used the Hybrid Monte Carlo updating algorithm. Each update consists 
in 10-15 leapfrog steps of At = 0.015 — 0.05 units of Langevin time followed 
by a Metropolis test with « 80% acceptance. We discarded approximately 
10000 initial sweeps for thermalization. The measures were separated by 
20-30 decorrelating sweeps. The errors were evaluated by the usual blocking 
procedure (see also later). 

We have generated 4 x 10 4 configurations of the reference field <p r in 
the thermalized regime. For a given (//"-configuration, by using the field 
transformation ([U]) from A[ to all considered Aj, we have computed the 
value of the remainder action 5Si and of 

& = 74 2 <l>n<f>m, 
^ nm 

and 

A*2i = - m) 2 <j) n <f) m . 

nm 

Finally, by taking the appropriate average, we obtain (see @) 

/x x \ Ej Xi exp{6Sj} gi exy{5Sj} 

X(A 2 ,A 4 )~— , // 2 (A 2 , A 4 ) ~ — . 19 

E 4 eocp{dSi} Ei exp{55i} 

for each A 2 , A4. Before illustrating the results, let us discuss the fluctuations 
of 5S. 



10 



3.3 Fluctuation of the remainder action 



We plot in Fig. ^] (the wider histograms) the distribution of the remainder 
action SSi for 8 4 , 16 4 and 20 4 lattices and for a typical value of A2 while 
A4 = A4. The fluctuations SSi are relatively small and therefore one expects 
that the statistical error of the original Monte Carlo does not substantially 
increase, even for the largest lattice. It is expected that the fluctuations grow 
as the square root of the lattice volume. Actually we find a slower growth 
up to L = 20. The width of the fluctuations for L = 8, L = 16 and L = 20 
are respectively 0.25, 0.75, 1.05 and therefore they appear to increase almost 
linearly with L. 

In view of possible future improvements it is interesting to examine which 
part of the remainder action gives the largest contribution to its fluctuations. 
We expand SS in powers of the field ( |T3D and analyze the effect of the first 
term. In Fig. we have plotted the distribution of the difference 

s's = ss- uc 2 - 1) £ 4>n ( Cf„ + C-„ ) • (20) 

nfi 

We see that the distributions in S'S are distinctly narrower than those of 
5S and so we conclude that, in this case, most of the SS fluctuations are 
due to its quadratic part. The importance of the quadratic terms in SS can 
(partially) be understood by observing that the field typically assumes small 
values. For the values of Aj considered for the distributions in Fig. |^, the 
values of <fi r n 2 and <\) 2 n are typically of the order of 0.4. Therefore the terms in 
the expansion (|14]) with higher powers in the fields are expected give only a 
small contribution, both to SS and to its fluctuations. 

For comparison, we have computed the fluctuations of the remainder 
action of the histogram method which is given by SS = J2 n [ v {4 > m K) ~ 
v((j) n , A[) ]. They are similar to those of the field transformation method as 
shown in Fig. |3|. Here we have plotted the widths of the distributions of the 
remainder action SS entering in the two methods together with the width of 
S'S in 
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Figure 2: Fluctuations of the remainder action for lattices of size L = 8, 16 
and 20. The field transformation leads from the reference point r (see Fig. 
to the point A2 = —0.2, A4 = 0.5. The wide histogram corresponds to the 



remainder action 5S (see (\ll)). The narrow histogram corresponds to 5'S, 
the remainder action minus the quadratic term (see ft2(\))- In this example 
the quartic coupling is the same, X2 = Xa- 
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Figure 3: The width of the fluctuations of the remainder action as a function 
of A2 for a L = 8 lattice. The reference point is r (see Fig. [![). The quartic 
coupling is A4 = 0.43 and thus A4 7^ A4. The open circles correspond to the 
fluctuations of5S (see fljj^)). The triangles stand for the fluctuations ofS'S 
(see The stars correspond to the fluctuations of AS = S[<f), A2, A4] — 

S[(j), Aj, A4] (the logarithm of the weight used in the histogram method). Errors 
have been estimated by naive blocking. 
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3.4 Discussion of errors 



The error analysis in the field transformation method is similar to that for the 
histogram method §, |. Let us assume that, after e.g. a standard blocking, 
the data are uncorrelated. First of all, one must take into consideration the 
correlation between numerator and denominator in ([19]). In general, for an 
observable T given by the ratio 



Ei fiWi _ (fw) 



(21) 



EiWi (g) 

with Wi a weight, the expression to be used for the propagation of errors is 



a 



fw 



(J , 



2a 2 t 

v fWW 



\ (fw) 2 (w) 2 (fw) (w) 



(22) 



where 



a 



°) w = 




2 _ 
°w = 


w - 


.2 _ 

fww 


(/» 2 > " 



(wY 



This expression explicitly shows that large fluctuations in the remainder ac- 
tion, giving large values of a w , would increase the statistical errors in the 
evaluation of the observables. This fact can be understood also by consider- 
ing that in this case only the few configurations corresponding to the large 
positive fluctuations of SS would contribute to the observables. 

A source of strong systematic error comes from the fact that we have 
sampled the Monte Carlo at the reference parameters A£ with a finite set of 
configurations. It is known |J that this error goes down with the logarithm 
of the number of configurations. 

Finally we recall that in this approach, measurements at different values 
of the parameters are obtained from the same set of configurations and there- 
fore all of them are strongly correlated. In practice the correlation matrix 
of the final measurements is almost singular. This problem is present in any 
kind of reweighting technique. 
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3.5 Results 



First we consider the 8 4 lattice. In Figs. we have plotted 1/%(A2, A4) at 
the values of parameters Aj corresponding to points in the three segments A, 
B and C (see Fig. [1]). The reported errors are obtained using (|22|). 

In order to check the accuracy of the results, we have performed, for 
some values of the parameters, independent Monte Carlo simulations with 
the same number of configurations. The corresponding results (black circles) 
are reported in Figs. ^-^| together with their statistical errors. The agreement 
between the results obtained by the field transformation and by the indepen- 
dent Monte Carlo runs is good even for values of the parameters not very 
close to the reference point r. We have also computed by the same technique 
the values of ^2(^2, A4) and found similar good agreement. 

In the various plots with fixed A4 (on the segments A, B and C) the 
values of 1/x decrease with A2, but they tend to flatten upward before the 
critical value is reached, as expected since we work with a finite lattice. We 
also find, as expected, that H2 is more sensitive to the finite volume effects 
(see (0)). 

The statistical errors obtained by the field transformation method are 
not really larger than those obtained by the direct Monte Carlo. This is due 
to the fact that, as shown in the previous subsection (see Figs. 0-0), the 
tfS'-fluctuations are not large. 

Results for larger volumes. Since the remainder action is a sum over the 
full lattice one may expect that the efficiency of the method will worsen by 
increasing the lattice volume. However we have found (see Figs. @-0) that, 
even for large lattices, the ^-fluctuations are not too large. Therefore the 
agreement with the results from a direct Monte Carlo test run should not be 
sizably degraded. 

This expectation is confirmed by our study for larger lattices. In Figs. |7| 
and H we have plotted the values of l/x(A 2 ,A 4 ) obtained for 16 4 and 20 4 
lattices respectively. Like in Figs. [||], we have reported also the measures 
obtained by Monte Carlo test simulations with the same number of indepen- 
dent configurations. The agreement with the direct Monte Carlo results is 
still good. As expected, the statistical errors are larger, since the fluctuations 
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Figure 4: Lattice of size L = 8. The emtpy circles represent the values of 
l/x(^2, A4) obtained by the field transformation method from the reference 
point r to values of (X2, ^a) in the segment A of Fig. H. The errors are purely 
statistical and they are computed according to §2^). We give also the values 
obtained by Monte Carlo test runs with the same statistics (black circles) and 
the values coming from the usual histogram method (stars). 
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Figure 5: Same as in Fig. f|, but for a mapping from the reference point r to 
values of (X2, X4) in the segment B of Fig. [J. 
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Figure 6: Same as in Fig. f|, but for a mapping from the reference point r 
values of (X2, X4) in the segment C of Fig. [J. 
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Figure 7: Lattice of size L—16. We have plotted l/x(A 2 , A 4 ) obtained by the 
field transformation method from the reference point r to parameters (X2, X4) 
in the segments B (see in Fig. [|). The errors are purely statistical and they 
are computed according to pi}/ We give also the values obtained by Monte 
Carlo test runs with the same statistics (black circles). 
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of the remainder action increase. However the increase of the errors is not 
dramatic and one still has a good determination of the observables. 

Comparison with the histogram method. In Figs. §-^we plot also the 
results obtained by the histogram method (stars). We have to distinguish 
the case in which the quartic coupling A4 is the same as that of the reference 
action A4 or is not. In the first case (see Fig. |5| in which A 4 = A4 = 0.5) 
we find that the results obtained by the Monte Carlo test simulations agree 
with those obtained both by the field transformation and by the histogram 
method. 

If instead A 4 7^ A4 (see Figs. £| and ||) we find that the results of the Monte 
Carlo test simulations disagree with the ones of the histogram method, while 
they still agree with the ones of the field transformation. This indicates that 
for A4 7^ A4 the regions of field configuration space probed by the histogram 
method for the two actions S[<p, AJ and S[4>, A£] have a small overlap. There- 
fore in order to obtain a reliable result by the histogram method one has to 
considerably increase the number of configurations sampled. The alternative 
we have proposed is to improve the overlapping by a field transformation. 

4 Discussion and conclusions 

The field transformation method, which we have introduced and analysed, 
can be viewed as a generalization of the usual histogram method. In principle 
it has the advantage that, by a suitable choice of the mapping, one has the 
possibility of reducing the remainder action, the crucial quantity entering into 
the weight which can endanger the statistical significance of the calculation. 
We have analyzed the case of a simple mapping, the ultralocal one in @. 
It gives a remainder action which has small fluctuations in a wide range of 
parameters and lattice sizes. Moreover the results for the susceptibility (and 
other computed quantities) obtained by the field transformation agree quite 
well with those obtained by Monte Carlo test runs, even when the lattice 
volume increases (at least up to L = 20, see Figs. £§-|]). This implies that 
the region of (//-configurations sampled by the Boltzmann weight e _s ^' ,A ^, 
when transformed by @, overlaps sufficiently well with the region probed 
by e~ s ^' Xi \ 
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Figure 8: Same as in Fig. Q/or a lattice of size L = 20. 
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Since the ultralocal transformation already gives small fluctuations and 
consistent results, we have not been forced to improve it. However we have 
indicated a possible strategy to further reduce the remainder action. The 
error analysis and the statistical independence of the results obtained by the 
field transformation method can be discussed analogously to the case of the 
histogram method 0, § . 

The histogram method corresponds to the trivial field transformation: 
<p r n = 4> n . For A4 7^ A4, we find that the results of this method disagree with 
those obtained by Monte Carlo test runs. This indicates that, for A 4 7^ A4, 
the relevant ^-configuration region for the action S[<f>, Aj] is not sufficiently 
probed by this method. 
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